This project seeks to estimate sport fish harvest and releases of rockfish in Alaska waters by improving on the Howard et al. (2020) methods and expand the time series back to 1977 when the statewide harvest survey (SWHS) was first implemented. This is essentially a Bayesian version of the Howard methods that allows for more appropriate and defensible sharing of information between areas, handles missing data in a more appropriate manor, accurately propagates uncertainty throughout the estimation procedure and replaces the Howard decision tree approach to low sample sizes with a hierarchical model. The methods and results for generating harvest estimates are generally consistent between the Bayesian model and the Howard methods. Harvest estimates are consistent with Howard estimates during contemporary times, but may differ based on more appropriate weighting of SWHS and logbook data, including estimating and correcting bias in the SWHS data.

Table 1. Summary of key improvements in reconstructiing sport fish removals of rockfish using the Bayesian model as compared to the Howard et al. (2020) methods.
Issue Howard Bayes
Time series 1999 - present 1977 - present
Bias in SWHS Not explicitly dealt with. Relies on logbook data and ratios of guided/unguided from SWHS data to estimate unguided releases and harvests. Explicitly estimates bias in SWHS harvest and release estimates based on logbook data.
Species composition of releases Assumes that species composition of releases is equal to that of the harvest, which is not evident in the logbook data. Recognizes different release probabilities by species / species assemblage and estimates it from logbook data and bias corrected SWHS data
Sample size limitations Uses sample size threshholds such that when areas fall below those threshholds values are borrowed from nearby areas. Uses a hierarchichacal modelling approach that shares information between areas in the same region. Thus all data is used, even with small sample sizes. This is a more sound method that avoids assumptions and uses all of the data.
Error propogation Error is propogated when variance estimates are available, but there is uncertainty associated with borrowing values from nearby areas, or the assumption of species compositions being identical in harvest and releases, are not dealt with. By breaking the assumption that species composition is equal between harvests and releases, uncertainty in the release estimates is more reflective of the fishery. Furthermore, the hyerarchichal approach more accurately captures uncertainy within and between areas within a region.
**Figure 1.**- Data sources for estimating rockfish harvests and releases in ADF&G commercial fisheries management units. Note that initial rockfish harvest estimates were not differentiated into species assemblage or species until 1998 when logbooks began differentiating by pelagic and non-pelagic. Logbooks began to collect data on yelloweye beginning in 2006. Port sampling programs to gather data on species composition of harvests began in 1996 in Southcentral and Kodiak and in 2006 in Southeast.

Figure 1.- Data sources for estimating rockfish harvests and releases in ADF&G commercial fisheries management units. Note that initial rockfish harvest estimates were not differentiated into species assemblage or species until 1998 when logbooks began differentiating by pelagic and non-pelagic. Logbooks began to collect data on yelloweye beginning in 2006. Port sampling programs to gather data on species composition of harvests began in 1996 in Southcentral and Kodiak and in 2006 in Southeast.


**Figure 2.**- SWHS harvest (left) and release (right) estimates from guided trips (x-axis) versus repoted harvests from charter logbooks (y-axis).

Figure 2.- SWHS harvest (left) and release (right) estimates from guided trips (x-axis) versus repoted harvests from charter logbooks (y-axis).


Results

**Figure X.**- Rhat values and proportion of parameters that converged (Rhat < 1.1.)

Figure X.- Rhat values and proportion of parameters that converged (Rhat < 1.1.)

Estimate comparison

Since previous estimates of rockfish harvest have been produced these first 3 graphs will be used to show how the modeled estimates compare to the estimates produced earlier. For total rockfish the estimates are in general agreement although differences are noted. These estimates should be more reliable because they include both SWHS and guide logbook data, handle variance more appropriately, use hierarchical distributions when data is missing, directly consider observation error and are produced using reproducible research.

**Figure 3.**- Total rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 3.- Total rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 3.**- Total rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 3.- Total rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.


Notes from Adam: When looking at only black rockfish the most significant differences are for the Prince William Sound Inside area. I did not spend a great deal of time tracking this down although it looks like the previous version used bad values for \(P_{(black)ayu}\) for at least unguided anglers. For the moment I would ignore the results for BSIA and SOKO2SAP. I think it is possible to give approximate values for these areas but it will require a little more coding which I have yet to do.

**Figure 4.**- Black rockfish harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 4.- Black rockfish harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.


And black rockfish releases…

**Figure 5.**- Black rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 5.- Black rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.





**Figure 6.**- Yellow rockfish harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 6.- Yellow rockfish harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 7.**- Yellow rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 7.- Yellow rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.





**Figure 8.**- DSR rockfish (excluding yelloweye) harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 8.- DSR rockfish (excluding yelloweye) harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 9.**- DSR rockfish releases (including yelloweye) 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 9.- DSR rockfish releases (including yelloweye) 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 11.**- Slope rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 11.- Slope rockfish harvests 1996-2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 12.**- Slope rockfish releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 12.- Slope rockfish releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

**Figure 13.**- *Other* rockfish (non-pelagic and non-yelloweye) harvests 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 13.- Other rockfish (non-pelagic and non-yelloweye) harvests 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



**Figure 14.**- *Other* rockfish (non-pelagic and non-yelloweye) releases 1996--2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.

Figure 14.- Other rockfish (non-pelagic and non-yelloweye) releases 1996–2023. Lines and error polygons represent model estimates and points and error bars represent Howard et al estimates.



Total Biomass Removal Estimates

**Figure 15.**- Black rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 15.- Black rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.



**Figure 16.**- Yellow rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 16.- Yellow rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

**Figure 17.**- Pelagic rockfish estimated total removals in lbs in 1996--2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 17.- Pelagic rockfish estimated total removals in lbs in 1996–2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


**Figure 18.**- Non-yelloweye, demersal shelf rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 18.- Non-yelloweye, demersal shelf rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


**Figure 19.**- Slope rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.

Figure 19.- Slope rockfish estimated total removals in lbs in 1996-2023. Total removals are calculated from the harvests, the estimated release mortality and mean fish weights. Error polygons represent 95% confidence intervals.


Model fit

Logbook residuals

**Figure 20.**- Residuals from logbook harvests.

Figure 20.- Residuals from logbook harvests.


SWHS residuals

**Figure 21.**- Residuals from SWHS harvests.

Figure 21.- Residuals from SWHS harvests.



**Figure 22.**- Residual of SWHS releases.

Figure 22.- Residual of SWHS releases.

Parameter estimates

P(Charter)

These histograms show the posterior distribution of the mean percent of rockfish harvested by the charter fleet.

**Figure 23.**- Mean percent of harvest by charter anglers.

Figure 23.- Mean percent of harvest by charter anglers.


When considered annually we see the percent of rockfish harvested by the charter fleet follows our data fairly well although the model smooths out the changes and we just do not have much information about this ratio. Prior to 2011 the percent charter is confounded with SWHS bias and should be mostly discounted.

**Figure 24.**- Annual estimates of the percent of harvest by charter anglers for 16 commerical fishing manamgent areas, 1996-2023.

Figure 24.- Annual estimates of the percent of harvest by charter anglers for 16 commerical fishing manamgent areas, 1996-2023.

P(Harvest)

These plots show the fitted logistic line to the proportion of caught rockfish that are harvested. These estimates are used for hindcasting catch estimates based on the harvest data in early years when catch estimates are unavailable.


**Figure 25.**- Annual proportion of pelagic rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.

Figure 25.- Annual proportion of pelagic rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.


**Figure 26.**- Annual proportion of yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.

Figure 26.- Annual proportion of yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available.


**Figure 27.**- Annual proportion of non-pelagic, non-yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note, that this is not estimated for Southeast areas because non=pelagics are divided between DSR (including yelloweye) and Slope species.

Figure 27.- Annual proportion of non-pelagic, non-yelloweye rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note, that this is not estimated for Southeast areas because non=pelagics are divided between DSR (including yelloweye) and Slope species.


**Figure 28.**- Annual proportion of DSR rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

Figure 28.- Annual proportion of DSR rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.


**Figure 29.**- Annual proportion of Slope rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

Figure 29.- Annual proportion of Slope rockfish catch that was harvested. Note that pre-1990 estimates are used to estimate catch in these years when catch estimates are not available. Note that the observed logbook data is for all non-pelagic, non-yelloweye fish.

SWHS bias

Figure 23 shows the mean estimate for SWHS bias in harvests and releases. Cook Inlet, North Gulf Coast and North Southeast Inside all look pretty good while most other areas have substantial bias. Prince William Sound Inside has the largest bias. Bias in release estimates is substantial and whereas the SWHS appears to underestimate harvests, it appears to greatly overestimates releases by a factor of 2 or more in most areas as derived from logbook reported releases.

**Figure 30.**- Mean SWHS bias for harvests and catches. Note that a bias < 1 indicates that the SWHS *underestimates* the true value and bias > 1 indicates the survey *overestimates* the true value.

Figure 30.- Mean SWHS bias for harvests and catches. Note that a bias < 1 indicates that the SWHS underestimates the true value and bias > 1 indicates the survey overestimates the true value.


Our estimates of SWHS harvest bias track observations fairly well when he have guided harvest estimates. The estimates of release bias in the SWHS data track observed patterns to an extent, but appear to smooth these more volatile disagreements with the logbook data. Adam postulated in his initial start on this that some of this could be the result of the estimates of the proportion guided. This value was not modelled with a trend and thus applies a constant estimate when hindcasting. Data on these relationships could greatly improve this model.

**Figure 31.**- Annual estimates of SWHS bias in harvests and releases for 16 commerical fishing manamgent areas, 1996-2023. Note that a bias < 1 indicates that the SWHS *underestimates* the true value and bias > 1 indicates the survey *overestimates* the true value.

Figure 31.- Annual estimates of SWHS bias in harvests and releases for 16 commerical fishing manamgent areas, 1996-2023. Note that a bias < 1 indicates that the SWHS underestimates the true value and bias > 1 indicates the survey overestimates the true value.

P(pelagic)

We model the percentage of pelagic rockfish in the harvest because we have the information for charter anglers (via logbooks) starting in 1998. Other than looking at the model estimates you can use Figure 25 to compare the two data streams for pelagic rockfish harvest. In general they are in agreement with major exceptions in Price William Sound inside, Prince William Sound outside (early in the time series) and South Southeast inside.

**Figure 32.**- Annual estimates of the percent of the sport harvest that was pelagic rockfish for 16 commerical fishing manamgent areas, 1996-2023.

Figure 32.- Annual estimates of the percent of the sport harvest that was pelagic rockfish for 16 commerical fishing manamgent areas, 1996-2023.

P(black|pelagic)

Note that in Southeast Alaska we only have composition data starting in 2006. Tania dug up old SE data, but it did not provide any useful data for species apportionment. For the most part, P(black|pelagic) is relatively constant across areas, with the exception of Cook Inlet and NSEI in Southeast AK. It may be worth discussing whether the shifts in those areas is a result of improved or changing species identification rather than actual shift in the species composition of the catch.

**Figure 33.**- Annual estimates of the percent of the sport harvest of pelagic rockfish that were black rockfish for 16 commerical fishing manamgent areas, 1996-2023. Kodiak panels include data from a hydroacoustic survey and the proportion of pelagic rockfish that are black in those areas (red) and the adjusted proportions based on obseved harvests for charter (blue) and private (cyan) users.

Figure 33.- Annual estimates of the percent of the sport harvest of pelagic rockfish that were black rockfish for 16 commerical fishing manamgent areas, 1996-2023. Kodiak panels include data from a hydroacoustic survey and the proportion of pelagic rockfish that are black in those areas (red) and the adjusted proportions based on obseved harvests for charter (blue) and private (cyan) users.

P(yelloweye|non-pelagic / yelloweye|DSR)

**Figure 34.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were yelloweye rockfish for 16 commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

Figure 34.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were yelloweye rockfish for 16 commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

P(DSR|non-pelagic)

**Figure 35.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were DSR rockfish for 6 Southeast commerical fishing manamgent areas, 1996-2023.

Figure 35.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were DSR rockfish for 6 Southeast commerical fishing manamgent areas, 1996-2023.

P(slope|non-pelagic)

**Figure 36.**- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were slope rockfish for 6 southeast commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.

Figure 36.- Annual estimates of the percent of the sport harvest of non-pelagic rockfish that were slope rockfish for 6 southeast commerical fishing manamgent areas, 1996-2023. Note that P(yelloweye) is the the proportion relative to non-pelagics for Central and Kodiak areas but is relative to DSR for Southeast areas.



P(slope|non-pelagic & non-yellowye) For release estimates

## NULL



Weight Fits

**Figure 37.**- Mean weights of black rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 37.- Mean weights of black rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 38.**- Mean weights of yelloweye rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 38.- Mean weights of yelloweye rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 39.**- Mean weights of non-black, pelagic rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 39.- Mean weights of non-black, pelagic rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 40.**- Mean weights of non-yelloweye, demersal shelf rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 40.- Mean weights of non-yelloweye, demersal shelf rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


**Figure 41.**- Mean weights of slope rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.

Figure 41.- Mean weights of slope rockfish used to estimate biomass from harvest and release mortality estimates in numbers of fish.


Summary of unconverged parameters:

The proportion of unconverged species composition proportional estimates that are near the bounds of 0 or 1 is 0.4150943

The proportion of unconverged true parameters (defined as Rhat > 1.01) that are also in excess of an Rhat of 1.1 is 0.08 while 0.92 are less than 1.1.

Table Xa. Species proportion derived quantities that are not with 0.05 of the 0 and 1 bounds and are not converged relative to an Rhat of 1.01.

variable

area_n

year

user

mean

median

lower

upper

Rhat

sd

skew

area

p_pelagic

11

4

2

0.397

0.392

0.182

0.646

1.011

0.121

0.041

CSEO

p_pelagic

11

7

1

0.465

0.464

0.225

0.711

1.013

0.126

0.008

CSEO

p_pelagic

13

1

1

0.255

0.233

0.077

0.546

1.079

0.123

0.179

NSEI

p_pelagic

13

1

2

0.100

0.083

0.026

0.260

1.076

0.063

0.271

NSEI

p_yellow

13

1

1

0.461

0.400

0.141

0.982

1.015

0.239

0.255

NSEI

p_yellow

13

1

2

0.463

0.401

0.139

0.983

1.018

0.238

0.260

NSEI

p_pelagic

13

10

1

0.281

0.264

0.095

0.562

1.068

0.122

0.140

NSEI

p_pelagic

13

10

2

0.114

0.100

0.032

0.282

1.061

0.065

0.216

NSEI

p_yellow

13

10

1

0.463

0.404

0.127

0.981

1.014

0.240

0.246

NSEI

p_yellow

13

10

2

0.463

0.403

0.137

0.984

1.013

0.241

0.249

NSEI

p_pelagic

13

11

1

0.289

0.273

0.104

0.551

1.045

0.120

0.133

NSEI

p_pelagic

13

11

2

0.115

0.100

0.030

0.285

1.058

0.067

0.223

NSEI

p_yellow

13

11

1

0.457

0.390

0.135

0.982

1.011

0.243

0.276

NSEI

p_yellow

13

11

2

0.464

0.406

0.130

0.984

1.012

0.241

0.240

NSEI

p_pelagic

13

12

1

0.289

0.272

0.101

0.559

1.044

0.121

0.140

NSEI

p_pelagic

13

12

2

0.118

0.103

0.034

0.287

1.062

0.067

0.225

NSEI

p_yellow

13

12

1

0.459

0.399

0.133

0.984

1.014

0.243

0.247

NSEI

p_yellow

13

12

2

0.466

0.405

0.135

0.981

1.013

0.239

0.255

NSEI

p_pelagic

13

13

1

0.299

0.283

0.106

0.581

1.058

0.123

0.130

NSEI

p_pelagic

13

13

2

0.122

0.109

0.035

0.294

1.045

0.068

0.192

NSEI

p_yellow

13

13

1

0.460

0.403

0.136

0.981

1.014

0.239

0.238

NSEI

p_yellow

13

13

2

0.465

0.410

0.136

0.983

1.014

0.239

0.230

NSEI

p_pelagic

13

14

1

0.284

0.270

0.108

0.553

1.038

0.113

0.124

NSEI

p_pelagic

13

14

2

0.121

0.110

0.036

0.279

1.048

0.063

0.175

NSEI

p_yellow

13

14

1

0.462

0.403

0.136

0.980

1.010

0.242

0.244

NSEI

p_yellow

13

14

2

0.464

0.404

0.133

0.981

1.018

0.240

0.250

NSEI

p_pelagic

13

15

1

0.260

0.241

0.096

0.524

1.034

0.110

0.173

NSEI

p_pelagic

13

15

2

0.118

0.105

0.035

0.272

1.041

0.063

0.207

NSEI

p_yellow

13

15

1

0.455

0.392

0.128

0.981

1.015

0.241

0.261

NSEI

p_yellow

13

15

2

0.458

0.395

0.131

0.983

1.013

0.242

0.260

NSEI

p_pelagic

13

16

1

0.300

0.288

0.115

0.558

1.033

0.116

0.104

NSEI

p_pelagic

13

16

2

0.130

0.118

0.040

0.287

1.029

0.064

0.186

NSEI

p_yellow

13

16

1

0.461

0.400

0.129

0.980

1.011

0.242

0.252

NSEI

p_yellow

13

16

2

0.461

0.401

0.133

0.982

1.016

0.240

0.250

NSEI

p_pelagic

13

17

1

0.316

0.305

0.127

0.574

1.030

0.115

0.096

NSEI

p_pelagic

13

17

2

0.139

0.127

0.042

0.311

1.031

0.069

0.175

NSEI

p_yellow

13

17

1

0.459

0.399

0.135

0.985

1.013

0.241

0.249

NSEI

p_yellow

13

17

2

0.463

0.405

0.139

0.983

1.016

0.238

0.244

NSEI

p_pelagic

13

18

1

0.325

0.317

0.139

0.579

1.019

0.114

0.070

NSEI

p_pelagic

13

18

2

0.142

0.130

0.047

0.305

1.017

0.068

0.177

NSEI

p_yellow

13

18

1

0.465

0.411

0.137

0.985

1.012

0.238

0.227

NSEI

p_yellow

13

18

2

0.461

0.402

0.137

0.982

1.016

0.238

0.248

NSEI

p_pelagic

13

19

1

0.328

0.318

0.137

0.578

1.017

0.115

0.087

NSEI

p_pelagic

13

19

2

0.147

0.134

0.050

0.306

1.023

0.068

0.191

NSEI

p_yellow

13

19

1

0.460

0.399

0.134

0.982

1.016

0.239

0.255

NSEI

p_yellow

13

19

2

0.459

0.396

0.130

0.983

1.012

0.242

0.260

NSEI

p_pelagic

13

2

1

0.261

0.238

0.081

0.552

1.058

0.124

0.186

NSEI

p_pelagic

13

2

2

0.102

0.086

0.027

0.266

1.094

0.063

0.255

NSEI

p_yellow

13

2

1

0.463

0.403

0.132

0.983

1.014

0.240

0.250

NSEI

p_yellow

13

2

2

0.464

0.406

0.142

0.981

1.012

0.239

0.243

NSEI

p_pelagic

13

20

1

0.346

0.333

0.156

0.587

1.018

0.114

0.114

NSEI

p_pelagic

13

20

2

0.155

0.142

0.054

0.332

1.011

0.073

0.178

NSEI

p_yellow

13

20

1

0.461

0.401

0.130

0.981

1.013

0.239

0.251

NSEI

p_yellow

13

20

2

0.466

0.408

0.140

0.983

1.012

0.238

0.244

NSEI

p_pelagic

13

21

1

0.365

0.355

0.170

0.610

1.012

0.117

0.085

NSEI

p_yellow

13

21

1

0.462

0.403

0.141

0.984

1.014

0.238

0.248

NSEI

p_yellow

13

21

2

0.463

0.406

0.132

0.982

1.013

0.238

0.239

NSEI

p_yellow

13

22

1

0.459

0.398

0.137

0.980

1.015

0.238

0.257

NSEI

p_yellow

13

22

2

0.462

0.402

0.138

0.983

1.015

0.236

0.254

NSEI

p_yellow

13

23

1

0.424

0.361

0.105

0.975

1.015

0.244

0.258

NSEI

p_yellow

13

23

2

0.508

0.467

0.162

0.986

1.017

0.230

0.178

NSEI

p_yellow

13

24

1

0.429

0.365

0.105

0.979

1.011

0.243

0.263

NSEI

p_yellow

13

25

1

0.436

0.378

0.110

0.978

1.017

0.244

0.238

NSEI

p_yellow

13

25

2

0.508

0.462

0.165

0.985

1.013

0.231

0.199

NSEI

p_yellow

13

26

1

0.422

0.362

0.108

0.977

1.013

0.238

0.252

NSEI

p_yellow

13

26

2

0.501

0.455

0.155

0.985

1.014

0.231

0.199

NSEI

p_yellow

13

27

1

0.409

0.346

0.100

0.977

1.013

0.238

0.265

NSEI

p_yellow

13

27

2

0.499

0.456

0.150

0.984

1.017

0.230

0.187

NSEI

p_yellow

13

28

1

0.424

0.368

0.113

0.970

1.013

0.230

0.243

NSEI

p_yellow

13

28

2

0.499

0.461

0.160

0.982

1.017

0.222

0.171

NSEI

p_yellow

13

29

1

0.390

0.337

0.110

0.934

1.018

0.216

0.245

NSEI

p_yellow

13

29

2

0.485

0.454

0.168

0.966

1.013

0.210

0.148

NSEI

p_pelagic

13

3

1

0.259

0.236

0.086

0.558

1.069

0.123

0.187

NSEI

p_pelagic

13

3

2

0.103

0.087

0.026

0.269

1.062

0.064

0.250

NSEI

p_yellow

13

3

1

0.460

0.399

0.139

0.980

1.013

0.236

0.258

NSEI

p_yellow

13

3

2

0.465

0.405

0.135

0.981

1.013

0.239

0.251

NSEI

p_pelagic

13

4

1

0.262

0.241

0.086

0.557

1.067

0.122

0.172

NSEI

p_pelagic

13

4

2

0.103

0.088

0.027

0.263

1.086

0.062

0.242

NSEI

p_yellow

13

4

1

0.462

0.401

0.131

0.979

1.016

0.238

0.257

NSEI

p_yellow

13

4

2

0.460

0.397

0.128

0.982

1.011

0.243

0.259

NSEI

p_pelagic

13

5

1

0.265

0.244

0.085

0.546

1.067

0.123

0.171

NSEI

p_pelagic

13

5

2

0.106

0.091

0.028

0.273

1.067

0.063

0.237

NSEI

p_yellow

13

5

1

0.464

0.403

0.139

0.982

1.014

0.238

0.257

NSEI

p_yellow

13

5

2

0.466

0.406

0.137

0.984

1.013

0.239

0.251

NSEI

p_pelagic

13

6

1

0.262

0.242

0.088

0.534

1.051

0.119

0.169

NSEI

p_pelagic

13

6

2

0.106

0.092

0.029

0.276

1.072

0.064

0.218

NSEI

p_yellow

13

6

1

0.462

0.406

0.132

0.981

1.013

0.239

0.235

NSEI

p_yellow

13

6

2

0.467

0.401

0.133

0.982

1.017

0.239

0.276

NSEI

p_pelagic

13

7

1

0.271

0.249

0.094

0.552

1.065

0.124

0.178

NSEI

p_pelagic

13

7

2

0.106

0.091

0.027

0.272

1.059

0.065

0.232

NSEI

p_yellow

13

7

1

0.458

0.401

0.129

0.981

1.013

0.241

0.237

NSEI

p_yellow

13

7

2

0.467

0.405

0.136

0.981

1.020

0.240

0.258

NSEI

p_pelagic

13

8

1

0.275

0.256

0.093

0.563

1.052

0.123

0.155

NSEI

p_pelagic

13

8

2

0.111

0.096

0.028

0.280

1.068

0.065

0.229

NSEI

p_yellow

13

8

1

0.463

0.403

0.131

0.980

1.013

0.240

0.250

NSEI

p_yellow

13

8

2

0.464

0.404

0.137

0.980

1.013

0.239

0.251

NSEI

p_pelagic

13

9

1

0.277

0.261

0.092

0.556

1.042

0.122

0.132

NSEI

p_pelagic

13

9

2

0.112

0.097

0.032

0.280

1.065

0.065

0.229

NSEI

p_yellow

13

9

1

0.462

0.402

0.134

0.983

1.012

0.239

0.251

NSEI

p_yellow

13

9

2

0.464

0.400

0.138

0.982

1.014

0.240

0.267

NSEI

p_yellow

14

11

1

0.578

0.564

0.255

0.984

1.014

0.185

0.076

NSEO

p_yellow

14

13

1

0.580

0.564

0.259

0.982

1.014

0.184

0.087

NSEO

p_yellow

14

13

2

0.620

0.613

0.299

0.988

1.014

0.177

0.040

NSEO

p_yellow

14

15

1

0.577

0.566

0.256

0.984

1.010

0.183

0.060

NSEO

p_yellow

14

15

2

0.620

0.611

0.291

0.988

1.010

0.177

0.051

NSEO

p_yellow

14

16

2

0.617

0.606

0.286

0.987

1.011

0.176

0.062

NSEO

p_yellow

14

19

1

0.583

0.568

0.261

0.984

1.010

0.184

0.082

NSEO

p_yellow

14

2

2

0.621

0.611

0.303

0.988

1.011

0.176

0.057

NSEO

p_yellow

14

20

1

0.582

0.567

0.268

0.985

1.011

0.182

0.082

NSEO

p_yellow

14

20

2

0.617

0.608

0.302

0.988

1.012

0.177

0.051

NSEO

p_yellow

14

21

2

0.614

0.606

0.291

0.985

1.011

0.177

0.045

NSEO

p_yellow

14

22

2

0.616

0.608

0.293

0.986

1.010

0.179

0.045

NSEO

p_yellow

14

23

2

0.612

0.604

0.282

0.988

1.014

0.182

0.044

NSEO

p_yellow

14

24

2

0.614

0.604

0.278

0.989

1.011

0.183

0.055

NSEO

p_yellow

14

25

2

0.616

0.610

0.278

0.988

1.011

0.179

0.033

NSEO

p_yellow

14

26

1

0.585

0.574

0.257

0.984

1.011

0.183

0.060

NSEO

p_yellow

14

3

2

0.616

0.605

0.292

0.988

1.011

0.180

0.061

NSEO

p_yellow

14

4

2

0.618

0.610

0.292

0.987

1.010

0.178

0.045

NSEO

p_yellow

14

5

2

0.617

0.614

0.284

0.988

1.010

0.179

0.017

NSEO

p_yellow

14

6

1

0.579

0.566

0.258

0.985

1.011

0.185

0.070

NSEO

p_black

14

7

1

0.918

0.945

0.690

0.993

1.014

0.084

-0.322

NSEO

p_yellow

14

8

1

0.577

0.562

0.252

0.983

1.010

0.186

0.081

NSEO

p_yellow

14

9

1

0.577

0.563

0.257

0.985

1.012

0.184

0.076

NSEO

p_yellow

14

9

2

0.618

0.611

0.290

0.987

1.010

0.179

0.039

NSEO

p_pelagic

16

1

1

0.393

0.386

0.143

0.677

1.014

0.142

0.049

SSEO

p_pelagic

16

10

1

0.395

0.387

0.158

0.671

1.013

0.137

0.058

SSEO

p_pelagic

16

10

2

0.375

0.368

0.141

0.646

1.018

0.134

0.052

SSEO

p_slope

16

10

2

0.076

0.057

0.012

0.261

1.013

0.065

0.290

SSEO

p_pelagic

16

11

1

0.395

0.391

0.156

0.675

1.011

0.136

0.029

SSEO

p_pelagic

16

11

2

0.375

0.368

0.144

0.641

1.010

0.132

0.053

SSEO

p_pelagic

16

15

2

0.379

0.372

0.148

0.652

1.010

0.131

0.054

SSEO

p_slope

16

17

2

0.078

0.059

0.010

0.265

1.017

0.067

0.283

SSEO

p_pelagic

16

2

2

0.371

0.366

0.132

0.654

1.012

0.137

0.036

SSEO

p_pelagic

16

3

1

0.391

0.383

0.143

0.668

1.016

0.139

0.057

SSEO

p_pelagic

16

3

2

0.371

0.366

0.135

0.649

1.013

0.136

0.037

SSEO

p_dsr

16

30

1

0.419

0.419

0.393

0.447

1.013

0.014

0.000

SSEO

p_dsr

16

31

1

0.401

0.400

0.374

0.432

1.098

0.015

0.067

SSEO

p_slope

16

36

1

0.061

0.061

0.041

0.079

1.026

0.009

0.000

SSEO

p_slope

16

37

1

0.065

0.066

0.047

0.082

1.011

0.009

-0.114

SSEO

p_pelagic

16

4

1

0.394

0.388

0.145

0.676

1.012

0.141

0.042

SSEO

p_pelagic

16

4

2

0.370

0.364

0.130

0.647

1.011

0.138

0.044

SSEO

p_pelagic

16

42

1

0.833

0.833

0.826

0.840

1.125

0.003

0.000

SSEO

p_yellow

16

42

1

0.406

0.407

0.385

0.429

1.044

0.011

-0.090

SSEO

p_dsr

16

42

1

0.451

0.442

0.404

0.509

3.457

0.030

0.298

SSEO

p_slope

16

42

1

0.167

0.180

0.110

0.209

4.883

0.032

-0.409

SSEO

p_pelagic

16

5

1

0.394

0.389

0.154

0.671

1.014

0.136

0.037

SSEO

p_pelagic

16

5

2

0.370

0.363

0.131

0.647

1.016

0.137

0.051

SSEO

p_pelagic

16

6

1

0.394

0.391

0.141

0.677

1.012

0.141

0.021

SSEO

p_pelagic

16

6

2

0.372

0.363

0.130

0.650

1.013

0.137

0.066

SSEO

p_pelagic

16

7

1

0.392

0.389

0.147

0.665

1.010

0.137

0.022

SSEO

p_pelagic

16

9

2

0.374

0.362

0.145

0.650

1.013

0.135

0.089

SSEO

p_pelagic

3

4

2

0.187

0.175

0.058

0.389

1.014

0.085

0.140

PWSI

p_pelagic

3

6

1

0.549

0.555

0.268

0.789

1.011

0.135

-0.045

PWSI

p_pelagic

3

7

2

0.188

0.176

0.062

0.388

1.011

0.086

0.140

PWSI

p_black

8

40

1

0.486

0.416

0.343

0.964

1.016

0.190

0.369

afognak

Table Xb. True parameters with an Rhat in excess of 1.1

variable

area_n

year/species

mean

median

lower

upper

Rhat

sd

skew

year

area

beta1_pH

11

3

82,647,383.439

6.859

0.216

694,725.175

1.188

2,584,012,155.858

0.032

CSEO

beta1_pH

12

3

605,770,439.507

6.632

0.227

2,316,460.836

1.276

24,991,659,588.711

0.024

EWYKT

beta1_pH

13

3

10,277,915,288.490

5.762

0.220

2,463,886.264

1.239

409,336,022,438.558

0.025

NSEI

beta0_pelagic

13

-1.285

-1.369

-1.965

-0.286

1.157

0.477

0.176

NSEI

beta1_pH

14

3

790,930,903.417

5.912

0.222

1,977,960.149

1.138

25,176,501,274.686

0.031

NSEO

beta1_pH

15

3

15,108,533,445.479

7.208

0.222

2,831,806.250

1.252

470,457,299,250.957

0.032

SSEI

beta1_pH

16

3

809,111,305,285.909

6.442

0.225

1,917,248.209

1.291

44,268,526,024,671.227

0.018

SSEO

beta3_pH

16

4

41.597

41.599

40.892

42.486

1.127

0.421

-0.005

SSEO

beta0_dsr

16

0.161

0.161

0.129

0.193

1.150

0.016

0.000

SSEO

tau_beta1_pH

3

4

0.758

0.087

0.011

4.752

1.250

3.293

0.204

PWSI

beta2_pH

9

3

0.190

0.160

0.082

0.458

1.149

0.224

0.134

eastside

Table Xc. Summary statistics describing the characteristics of the derived harvest and release estimates that are not converged relative to an Rhat of 1.01.

Statistic

non.BSAI.SOKO.Harv

BSAI.SOKO.Harv

non.BSAI.SOKO.Rel

BSAI.SOKO.Rel

n unconverged at Rhat = 1.01

278.00

214.00

825.00

330.00

Prop. of total

0.57

0.43

0.71

0.29

Prop. median less than 1

0.36

0.26

0.35

0.42

Prop. median less than 10

0.46

0.60

0.56

0.76

Prop. median less than 50

0.52

0.79

0.70

0.84

Median of the median values

33.89

5.04

4.66

1.77

Mean of the median values

684.43

250.44

273.32

82.56

Mean skew of posteriors

0.24

0.28

0.27

0.24

Table Xd. Summary statistics describing the characteristics of the derived harvest and release estimates that are not converged relative to an Rhat of 1.1.

Statistic

non.BSAI.SOKO.Harv

BSAI.SOKO.Harv

non.BSAI.SOKO.Rel

BSAI.SOKO.Rel

n unconverged at Rhat = 1.1

5.00

8.00

39.00

34.00

Prop. of total

0.38

0.62

0.53

0.47

Prop. median less than 1

0.80

0.50

0.90

0.56

Prop. median less than 10

1.00

0.88

0.95

0.94

Prop. median less than 50

1.00

1.00

0.95

0.94

Median of the median values

0.00

1.96

0.01

0.48

Mean of the median values

1.50

3.88

39.18

20.57

Mean skew of posteriors

0.17

0.21

0.15

0.17

Table Xe. All parameters and derived quantities that are not converged relative to an Rhat of 1.1.

variable

Freq

beta0_dsr

1

beta0_pelagic

1

beta1_pH

6

beta2_pH

1

beta3_pH

1

eps_dsr

1

eps_pH

6

eps_sl

1

H_ayg

1

H_ayu

1

Hb_ayg

1

Hb_ayu

1

Hd_ay

1

Hd_ayg

2

Ho_ayu

5

Hp_ayg

1

Hp_ayu

1

Hs_ay

2

Hs_ayg

2

Hs_ayu

7

Hy_ay

1

Hy_ayg

2

Hy_ayu

1

p_black

1

p_dsr

1

p_pelagic

1

p_slope

2

pH

14

R_ay

1

R_ayg

1

R_ayu

2

Rb_ay

1

Rb_ayg

1

Rb_ayu

1

Rd_ay

6

Rd_ayg

5

Rd_ayu

4

re_dsr

2

re_pH

9

re_slope

2

Ro_ay

13

Ro_ayg

6

Ro_ayu

32

Rp_ay

1

Rp_ayg

1

Rp_ayu

1

Rs_ay

7

Rs_ayg

13

Rs_ayu

22

Ry_ay

5

Ry_ayg

3

Ry_ayu

15

tau_beta1_pH

1

Parameter estimates:

Summary Table of Parameter Estimates
Parameter mean sd Lower_CI Median Upper_CI
mu_bc_H[1] -0.016 0.071 -0.148 -0.019 0.134
mu_bc_H[2] 0.075 0.054 -0.024 0.072 0.187
mu_bc_H[3] -0.273 0.067 -0.401 -0.274 -0.142
mu_bc_H[4] -0.918 0.180 -1.272 -0.916 -0.567
mu_bc_H[5] 1.337 1.169 0.097 1.061 4.347
mu_bc_H[6] -2.364 0.316 -2.957 -2.376 -1.724
mu_bc_H[7] -0.452 0.102 -0.661 -0.448 -0.258
mu_bc_H[8] 0.265 0.351 -0.306 0.231 1.075
mu_bc_H[9] -0.245 0.125 -0.481 -0.248 0.010
mu_bc_H[10] 0.060 0.064 -0.055 0.058 0.191
mu_bc_H[11] -0.125 0.037 -0.199 -0.125 -0.053
mu_bc_H[12] -0.249 0.106 -0.483 -0.241 -0.051
mu_bc_H[13] -0.163 0.072 -0.305 -0.163 -0.020
mu_bc_H[14] -0.306 0.092 -0.503 -0.302 -0.137
mu_bc_H[15] -0.339 0.047 -0.428 -0.341 -0.243
mu_bc_H[16] -0.156 0.368 -0.784 -0.181 0.655
mu_bc_R[1] 1.042 0.121 0.803 1.044 1.281
mu_bc_R[2] 1.015 0.076 0.857 1.017 1.157
mu_bc_R[3] 0.219 0.091 0.039 0.221 0.393
mu_bc_R[4] 0.666 0.197 0.267 0.674 1.030
mu_bc_R[5] 1.673 0.368 0.938 1.672 2.384
mu_bc_R[6] -1.410 0.345 -2.092 -1.408 -0.754
mu_bc_R[7] 0.515 0.169 0.176 0.525 0.827
mu_bc_R[8] 0.643 0.147 0.357 0.646 0.932
mu_bc_R[9] 0.379 0.174 0.011 0.385 0.701
mu_bc_R[10] 0.895 0.104 0.687 0.895 1.096
mu_bc_R[11] 1.108 0.076 0.957 1.107 1.259
mu_bc_R[12] 0.872 0.165 0.546 0.874 1.189
mu_bc_R[13] 1.041 0.091 0.858 1.041 1.218
mu_bc_R[14] 0.860 0.122 0.615 0.860 1.099
mu_bc_R[15] 0.796 0.097 0.601 0.796 0.984
mu_bc_R[16] 1.347 0.113 1.125 1.349 1.564
tau_pH[1] 2.448 0.183 2.108 2.443 2.827
tau_pH[2] 1.571 0.147 1.295 1.568 1.875
tau_pH[3] 2.273 0.300 1.739 2.256 2.922
tau_pH[4] 0.991 0.137 0.751 0.982 1.282
tau_pH[5] 0.425 0.075 0.301 0.419 0.592
beta0_pH[1,1] 0.548 0.176 0.181 0.558 0.865
beta0_pH[2,1] 1.403 0.183 0.988 1.416 1.717
beta0_pH[3,1] 1.390 0.232 0.821 1.420 1.760
beta0_pH[4,1] 1.541 0.280 0.841 1.578 1.962
beta0_pH[5,1] -0.814 0.398 -1.810 -0.771 -0.152
beta0_pH[6,1] -0.607 0.522 -1.995 -0.456 0.038
beta0_pH[7,1] 0.752 0.011 0.732 0.752 0.773
beta0_pH[8,1] -0.699 0.411 -1.809 -0.643 -0.063
beta0_pH[9,1] -0.343 0.297 -0.990 -0.313 0.190
beta0_pH[10,1] 0.262 0.206 -0.179 0.276 0.631
beta0_pH[11,1] 0.076 0.274 -0.488 0.075 0.614
beta0_pH[12,1] 0.614 0.195 0.225 0.618 0.991
beta0_pH[13,1] -0.651 0.154 -0.970 -0.645 -0.367
beta0_pH[14,1] 0.010 0.178 -0.348 0.011 0.350
beta0_pH[15,1] 0.050 0.203 -0.443 0.069 0.372
beta0_pH[16,1] -0.843 0.205 -1.262 -0.836 -0.458
beta0_pH[1,2] 2.566 0.084 2.403 2.565 2.732
beta0_pH[2,2] 3.079 0.022 3.035 3.079 3.123
beta0_pH[3,2] 3.626 0.041 3.548 3.625 3.704
beta0_pH[4,2] 2.936 0.019 2.899 2.936 2.974
beta0_pH[5,2] 4.603 1.531 2.762 4.261 8.611
beta0_pH[6,2] 3.331 0.161 3.027 3.330 3.657
beta0_pH[7,2] 2.026 0.049 1.930 2.025 2.120
beta0_pH[8,2] 2.822 0.065 2.696 2.821 2.950
beta0_pH[9,2] 3.962 0.210 3.563 3.959 4.387
beta0_pH[10,2] 4.002 0.129 3.757 4.002 4.266
beta0_pH[11,2] 1.885 0.009 1.868 1.885 1.903
beta0_pH[12,2] 1.801 0.069 1.669 1.801 1.937
beta0_pH[13,2] 2.497 0.027 2.444 2.497 2.548
beta0_pH[14,2] 1.803 0.024 1.758 1.804 1.850
beta0_pH[15,2] 2.545 0.017 2.512 2.545 2.577
beta0_pH[16,2] 2.765 0.025 2.716 2.765 2.814
beta0_pH[1,3] 0.698 0.581 -0.869 0.854 1.378
beta0_pH[2,3] 2.126 0.017 2.092 2.126 2.160
beta0_pH[3,3] 2.689 0.024 2.642 2.689 2.738
beta0_pH[4,3] 2.981 0.035 2.913 2.980 3.047
beta0_pH[5,3] 3.769 2.556 0.977 3.110 10.749
beta0_pH[6,3] 1.258 0.325 0.701 1.238 1.953
beta0_pH[7,3] 0.857 0.050 0.759 0.858 0.958
beta0_pH[8,3] 0.285 0.124 0.028 0.293 0.509
beta0_pH[9,3] -2.816 1.047 -5.002 -2.726 -0.975
beta0_pH[10,3] -0.885 0.908 -2.981 -0.758 0.568
beta0_pH[11,3] 0.199 11.452 -22.417 0.169 22.912
beta0_pH[12,3] 0.171 11.285 -21.860 0.145 22.611
beta0_pH[13,3] 0.121 11.558 -21.786 0.088 23.416
beta0_pH[14,3] 0.131 11.526 -22.439 0.143 22.457
beta0_pH[15,3] -0.020 11.428 -22.742 0.108 21.327
beta0_pH[16,3] 0.093 11.312 -21.325 0.077 21.749
beta0_pH[11,4] 1.104 0.464 0.200 1.117 1.972
beta0_pH[12,4] 0.953 0.735 -0.649 1.005 2.308
beta0_pH[13,4] 1.051 0.565 -0.315 1.147 1.761
beta0_pH[14,4] 0.507 0.580 -0.764 0.559 1.475
beta0_pH[15,4] 0.993 0.461 -0.009 1.062 1.682
beta0_pH[16,4] 0.899 0.605 -0.416 0.940 1.974
beta0_pH[11,5] 2.973 0.147 2.698 2.970 3.270
beta0_pH[12,5] -1.026 0.456 -1.912 -1.017 -0.135
beta0_pH[13,5] 4.914 0.412 4.138 4.897 5.733
beta0_pH[14,5] 3.904 0.358 3.239 3.891 4.627
beta0_pH[15,5] 2.292 0.124 2.057 2.289 2.534
beta0_pH[16,5] 0.250 0.063 0.160 0.238 0.400
beta1_pH[1,1] 2.591 0.490 1.816 2.529 3.810
beta1_pH[2,1] 1.676 0.438 1.062 1.608 2.888
beta1_pH[3,1] 1.824 0.470 1.113 1.756 3.115
beta1_pH[4,1] 1.998 0.535 1.282 1.880 3.544
beta1_pH[5,1] 2.347 0.531 1.519 2.276 3.617
beta1_pH[6,1] 2.296 0.915 1.097 2.040 4.572
beta1_pH[7,1] 0.000 0.000 0.000 0.000 0.000
beta1_pH[8,1] 2.645 1.039 1.269 2.340 4.815
beta1_pH[9,1] 2.142 0.684 1.251 2.003 4.176
beta1_pH[10,1] 2.176 0.409 1.530 2.130 3.095
beta1_pH[11,1] 3.315 0.459 2.476 3.280 4.304
beta1_pH[12,1] 2.919 0.287 2.362 2.915 3.492
beta1_pH[13,1] 3.444 0.245 2.994 3.432 3.946
beta1_pH[14,1] 2.853 0.329 2.246 2.847 3.555
beta1_pH[15,1] 2.408 0.291 1.905 2.384 3.102
beta1_pH[16,1] 3.989 0.288 3.433 3.984 4.570
beta1_pH[1,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[2,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[3,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[4,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[5,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[6,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[7,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[8,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[9,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[10,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[11,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[12,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[13,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[14,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[15,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[16,2] 0.000 0.000 0.000 0.000 0.000
beta1_pH[1,3] 3.421 1.667 1.217 2.972 7.420
beta1_pH[2,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[3,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[4,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[5,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[6,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[7,3] 0.000 0.000 0.000 0.000 0.000
beta1_pH[8,3] 4.973 2.350 2.180 4.204 10.991
beta1_pH[9,3] 8.336 2.546 4.135 7.992 13.841
beta1_pH[10,3] 7.198 2.746 3.248 6.670 13.530
beta1_pH[11,3] 82647383.439 2584012155.858 0.216 6.859 694725.175
beta1_pH[12,3] 605770439.507 24991659588.711 0.227 6.632 2316460.836
beta1_pH[13,3] 10277915288.490 409336022438.558 0.220 5.762 2463886.264
beta1_pH[14,3] 790930903.417 25176501274.686 0.222 5.912 1977960.149
beta1_pH[15,3] 15108533445.479 470457299250.957 0.222 7.208 2831806.250
beta1_pH[16,3] 809111305285.909 44268526024671.227 0.225 6.442 1917248.209
beta1_pH[11,4] 1.665 0.550 0.596 1.663 2.705
beta1_pH[12,4] 3.191 0.745 1.802 3.143 4.834
beta1_pH[13,4] 0.910 0.595 0.130 0.828 2.297
beta1_pH[14,4] 1.836 0.669 0.646 1.792 3.256
beta1_pH[15,4] 0.954 0.532 0.137 0.900 2.116
beta1_pH[16,4] 3.546 0.689 2.308 3.513 5.042
beta1_pH[11,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[12,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[13,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[14,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[15,5] 0.000 0.000 0.000 0.000 0.000
beta1_pH[16,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[1,1] 1.207 0.850 0.323 0.897 3.343
beta2_pH[2,1] 1.540 0.951 0.188 1.345 3.402
beta2_pH[3,1] 1.347 0.842 0.188 1.187 2.925
beta2_pH[4,1] 1.051 0.791 0.169 0.769 2.881
beta2_pH[5,1] 2.135 1.800 0.255 1.654 7.177
beta2_pH[6,1] 0.944 1.482 0.104 0.371 5.684
beta2_pH[7,1] 0.000 0.000 0.000 0.000 0.000
beta2_pH[8,1] 0.717 1.185 0.108 0.313 4.450
beta2_pH[9,1] 0.972 1.356 0.143 0.494 5.479
beta2_pH[10,1] 1.251 1.304 0.254 0.812 5.117
beta2_pH[11,1] 0.429 0.240 0.203 0.367 1.028
beta2_pH[12,1] 1.537 0.786 0.615 1.296 3.586
beta2_pH[13,1] 1.016 0.417 0.533 0.912 2.146
beta2_pH[14,1] 1.123 0.674 0.421 0.910 3.266
beta2_pH[15,1] 1.034 0.617 0.337 0.863 2.957
beta2_pH[16,1] 1.063 0.376 0.634 0.983 1.952
beta2_pH[1,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[2,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[3,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[4,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[5,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[6,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[7,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[8,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[9,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[10,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[11,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[12,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[13,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[14,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[15,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[16,2] 0.000 0.000 0.000 0.000 0.000
beta2_pH[1,3] 2.381 4.195 0.114 0.517 15.486
beta2_pH[2,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[3,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[4,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[5,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[6,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[7,3] 0.000 0.000 0.000 0.000 0.000
beta2_pH[8,3] 2.226 1.797 0.584 1.481 7.092
beta2_pH[9,3] 0.190 0.224 0.082 0.160 0.458
beta2_pH[10,3] 0.302 0.358 0.098 0.231 0.860
beta2_pH[11,3] -8.416 6.052 -22.751 -7.268 -0.361
beta2_pH[12,3] -8.320 5.983 -22.509 -7.167 -0.387
beta2_pH[13,3] -8.238 6.015 -23.060 -6.960 -0.379
beta2_pH[14,3] -8.300 6.057 -23.037 -7.043 -0.445
beta2_pH[15,3] -8.258 5.985 -22.715 -6.984 -0.360
beta2_pH[16,3] -8.271 6.059 -22.872 -7.028 -0.366
beta2_pH[11,4] -7.325 4.922 -18.726 -6.270 -0.975
beta2_pH[12,4] -7.608 4.879 -18.840 -6.613 -1.210
beta2_pH[13,4] -7.351 4.895 -18.960 -6.262 -0.997
beta2_pH[14,4] -7.288 4.890 -18.740 -6.112 -1.134
beta2_pH[15,4] -7.375 4.835 -18.653 -6.312 -1.025
beta2_pH[16,4] -7.493 4.741 -18.623 -6.330 -1.492
beta2_pH[11,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[12,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[13,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[14,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[15,5] 0.000 0.000 0.000 0.000 0.000
beta2_pH[16,5] 0.000 0.000 0.000 0.000 0.000
beta3_pH[1,1] 34.968 1.170 32.859 34.955 37.465
beta3_pH[2,1] 32.957 1.554 30.448 32.786 36.768
beta3_pH[3,1] 33.052 1.542 30.185 32.978 36.858
beta3_pH[4,1] 32.816 2.058 29.480 32.644 37.828
beta3_pH[5,1] 27.638 1.291 25.958 27.479 30.657
beta3_pH[6,1] 33.160 3.434 26.030 33.071 41.094
beta3_pH[7,1] 25.000 0.000 25.000 25.000 25.000
beta3_pH[8,1] 35.525 3.952 28.812 35.581 42.722
beta3_pH[9,1] 32.459 3.160 27.678 31.899 40.561
beta3_pH[10,1] 32.583 1.355 30.164 32.498 35.397
beta3_pH[11,1] 31.264 1.109 29.311 31.210 33.534
beta3_pH[12,1] 30.763 0.480 29.770 30.795 31.637
beta3_pH[13,1] 31.414 0.547 30.397 31.398 32.588
beta3_pH[14,1] 32.523 0.770 31.248 32.463 34.160
beta3_pH[15,1] 30.953 0.759 29.372 30.970 32.423
beta3_pH[16,1] 30.147 0.356 29.434 30.148 30.826
beta3_pH[1,2] 34.485 5.403 25.545 34.299 43.371
beta3_pH[2,2] 34.471 5.489 25.438 34.450 43.534
beta3_pH[3,2] 30.332 7.227 18.606 29.980 42.343
beta3_pH[4,2] 30.480 7.103 18.683 30.478 42.397
beta3_pH[5,2] 30.454 7.235 18.692 30.296 42.381
beta3_pH[6,2] 30.286 7.178 18.556 30.391 42.300
beta3_pH[7,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[8,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[9,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[10,2] 25.000 0.000 25.000 25.000 25.000
beta3_pH[11,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[12,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[13,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[14,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[15,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[16,2] 24.000 0.000 24.000 24.000 24.000
beta3_pH[1,3] 39.767 2.580 33.478 40.165 43.704
beta3_pH[2,3] 34.590 5.457 25.497 34.666 43.536
beta3_pH[3,3] 30.342 7.185 18.719 30.333 42.369
beta3_pH[4,3] 30.434 7.208 18.572 30.474 42.391
beta3_pH[5,3] 25.000 0.000 25.000 25.000 25.000
beta3_pH[6,3] 25.000 0.000 25.000 25.000 25.000
beta3_pH[7,3] 25.000 0.000 25.000 25.000 25.000
beta3_pH[8,3] 42.178 0.842 40.971 41.918 43.854
beta3_pH[9,3] 35.943 1.690 34.522 35.306 40.930
beta3_pH[10,3] 37.788 2.468 34.615 37.174 43.333
beta3_pH[11,3] 38.440 2.598 34.205 38.467 42.768
beta3_pH[12,3] 38.469 2.585 34.261 38.452 42.777
beta3_pH[13,3] 38.506 2.619 34.298 38.500 42.782
beta3_pH[14,3] 38.493 2.575 34.218 38.512 42.766
beta3_pH[15,3] 38.459 2.579 34.255 38.402 42.766
beta3_pH[16,3] 38.499 2.577 34.169 38.523 42.757
beta3_pH[11,4] 40.358 1.095 38.078 40.318 42.690
beta3_pH[12,4] 41.877 0.827 40.224 42.186 42.908
beta3_pH[13,4] 41.563 1.463 38.029 41.737 43.843
beta3_pH[14,4] 40.387 0.846 38.583 40.427 41.881
beta3_pH[15,4] 40.988 1.381 37.752 40.977 43.611
beta3_pH[16,4] 41.597 0.421 40.892 41.599 42.486
beta3_pH[11,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[12,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[13,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[14,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[15,5] 25.000 0.000 25.000 25.000 25.000
beta3_pH[16,5] 25.000 0.000 25.000 25.000 25.000
beta4_pH[1,1] 0.262 0.021 0.223 0.262 0.302
beta4_pH[2,1] 0.211 0.005 0.201 0.211 0.221
beta4_pH[3,1] 0.173 0.007 0.158 0.173 0.186
beta4_pH[4,1] 0.294 0.011 0.271 0.294 0.317
beta4_pH[5,1] 0.415 0.253 0.082 0.376 0.974
beta4_pH[6,1] 0.767 0.375 0.322 0.728 1.448
beta4_pH[7,1] 0.849 0.346 0.353 0.795 1.647
beta4_pH[8,1] 0.155 0.069 0.042 0.148 0.303
beta4_pH[9,1] 0.214 0.088 0.059 0.210 0.402
beta4_pH[10,1] 0.001 0.001 0.000 0.001 0.003
beta4_pH[11,1] 0.446 0.012 0.424 0.446 0.470
beta4_pH[12,1] 0.642 0.022 0.600 0.641 0.686
beta4_pH[13,1] 0.631 0.013 0.606 0.631 0.658
beta4_pH[14,1] 0.669 0.031 0.611 0.669 0.735
beta4_pH[15,1] 0.797 0.021 0.757 0.797 0.840
beta4_pH[16,1] 1.223 0.041 1.147 1.222 1.304
beta4_pH[1,2] 0.294 0.096 0.120 0.291 0.499
beta4_pH[2,2] 0.514 0.023 0.471 0.515 0.558
beta4_pH[3,2] 0.446 0.012 0.423 0.445 0.469
beta4_pH[4,2] 0.654 0.037 0.581 0.653 0.728
beta4_pH[5,2] 0.723 0.838 0.113 0.533 2.595
beta4_pH[6,2] 0.903 0.764 0.206 0.732 2.664
beta4_pH[7,2] 0.885 0.421 0.357 0.813 1.886
beta4_pH[8,2] 0.853 0.207 0.474 0.850 1.281
beta4_pH[9,2] 1.152 1.146 0.295 0.870 3.852
beta4_pH[10,2] 0.457 0.052 0.358 0.456 0.563
beta4_pH[11,2] 1.269 0.050 1.172 1.269 1.371
beta4_pH[12,2] 0.604 0.111 0.401 0.601 0.832
beta4_pH[13,2] 0.946 0.045 0.861 0.945 1.037
beta4_pH[14,2] 0.935 0.062 0.820 0.932 1.060
beta4_pH[15,2] 1.012 0.025 0.965 1.011 1.063
beta4_pH[16,2] 1.409 0.076 1.262 1.410 1.560
beta4_pH[1,3] 0.240 0.076 0.102 0.236 0.397
beta4_pH[2,3] 0.050 0.017 0.020 0.049 0.085
beta4_pH[3,3] 0.052 0.010 0.033 0.052 0.072
beta4_pH[4,3] 0.167 0.032 0.104 0.167 0.229
beta4_pH[5,3] 0.779 1.006 0.111 0.565 2.837
beta4_pH[6,3] 0.894 0.724 0.198 0.741 2.702
beta4_pH[7,3] 0.887 0.450 0.355 0.815 1.814
beta4_pH[8,3] 0.393 0.188 0.106 0.368 0.833
beta4_pH[9,3] 0.924 0.866 0.162 0.706 2.776
beta4_pH[10,3] 0.056 0.038 0.003 0.050 0.146
beta4_pH[11,3] 0.637 0.233 0.297 0.600 1.222
beta4_pH[12,3] 0.692 0.219 0.362 0.658 1.219
beta4_pH[13,3] 0.706 0.177 0.426 0.682 1.115
beta4_pH[14,3] 0.870 0.220 0.520 0.842 1.374
beta4_pH[15,3] 0.816 0.225 0.468 0.789 1.335
beta4_pH[16,3] 1.391 1.745 0.122 0.873 5.542
beta4_pH[1,4] 0.632 0.381 0.195 0.555 1.540
beta4_pH[2,4] 0.804 0.890 0.118 0.581 2.742
beta4_pH[3,4] 0.872 0.921 0.116 0.598 3.378
beta4_pH[4,4] 0.746 0.569 0.203 0.617 2.044
beta4_pH[5,4] 0.781 0.887 0.114 0.569 2.746
beta4_pH[6,4] 0.911 0.789 0.204 0.729 2.720
beta4_pH[7,4] 0.886 0.417 0.353 0.807 1.972
beta4_pH[8,4] 0.731 0.631 0.151 0.578 2.056
beta4_pH[9,4] 0.745 0.889 0.111 0.552 2.679
beta4_pH[10,4] 1.424 2.960 0.007 0.373 11.185
beta4_pH[11,4] 0.445 0.017 0.412 0.445 0.479
beta4_pH[12,4] 0.562 0.053 0.464 0.560 0.668
beta4_pH[13,4] 0.925 0.036 0.853 0.924 0.995
beta4_pH[14,4] 0.922 0.074 0.782 0.919 1.072
beta4_pH[15,4] 0.715 0.020 0.676 0.715 0.754
beta4_pH[16,4] 0.686 0.042 0.607 0.684 0.773
beta4_pH[1,5] 0.625 0.368 0.199 0.542 1.516
beta4_pH[2,5] 0.766 0.729 0.106 0.564 2.742
beta4_pH[3,5] 0.845 0.877 0.120 0.604 3.153
beta4_pH[4,5] 0.725 0.481 0.197 0.617 1.912
beta4_pH[5,5] 0.763 0.947 0.114 0.562 2.680
beta4_pH[6,5] 0.903 0.828 0.204 0.723 2.652
beta4_pH[7,5] 0.884 0.432 0.345 0.800 1.941
beta4_pH[8,5] 0.737 0.626 0.145 0.586 2.119
beta4_pH[9,5] 0.790 1.038 0.120 0.556 2.852
beta4_pH[10,5] 1.521 3.055 0.006 0.377 11.118
beta4_pH[11,5] 0.606 0.053 0.504 0.606 0.715
beta4_pH[12,5] 0.617 0.170 0.332 0.601 0.990
beta4_pH[13,5] 0.772 0.061 0.664 0.769 0.899
beta4_pH[14,5] 0.847 0.123 0.629 0.841 1.104
beta4_pH[15,5] 1.142 0.075 1.000 1.141 1.299
beta4_pH[16,5] 18.829 3.879 11.092 19.175 24.707
beta0_pelagic[1] 2.268 0.036 2.201 2.267 2.342
beta0_pelagic[2] 1.537 0.007 1.524 1.538 1.551
beta0_pelagic[3] 0.173 0.277 -0.545 0.239 0.557
beta0_pelagic[4] 0.241 0.197 -0.186 0.252 0.612
beta0_pelagic[5] 1.449 0.159 1.142 1.447 1.768
beta0_pelagic[6] 1.355 0.061 1.243 1.355 1.478
beta0_pelagic[7] 1.471 0.038 1.399 1.470 1.550
beta0_pelagic[8] 1.672 0.043 1.591 1.671 1.757
beta0_pelagic[9] 2.774 0.062 2.659 2.772 2.900
beta0_pelagic[10] 2.603 0.032 2.544 2.602 2.667
beta0_pelagic[11] -0.140 0.169 -0.541 -0.127 0.149
beta0_pelagic[12] 1.800 0.035 1.734 1.799 1.870
beta0_pelagic[13] -1.285 0.477 -1.965 -1.369 -0.286
beta0_pelagic[14] -0.284 0.187 -0.706 -0.260 0.016
beta0_pelagic[15] -0.209 0.035 -0.303 -0.195 -0.180
beta0_pelagic[16] -0.496 0.405 -1.392 -0.401 0.068
beta1_pelagic[1] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[2] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[3] 0.782 0.311 0.318 0.722 1.643
beta1_pelagic[4] 0.804 0.221 0.387 0.797 1.258
beta1_pelagic[5] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[6] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[7] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[8] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[9] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[10] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[11] 1.148 0.239 0.736 1.126 1.701
beta1_pelagic[12] 0.000 0.000 0.000 0.000 0.000
beta1_pelagic[13] 3.356 0.771 1.806 3.286 4.838
beta1_pelagic[14] 1.771 0.432 1.110 1.690 2.816
beta1_pelagic[15] 0.082 0.150 0.000 0.012 0.577
beta1_pelagic[16] 1.658 0.633 0.759 1.510 2.970
beta2_pelagic[1] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[2] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[3] 6.125 5.970 0.180 4.507 21.249
beta2_pelagic[4] 7.492 5.272 0.548 6.607 20.718
beta2_pelagic[5] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[6] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[7] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[8] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[9] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[10] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[11] 1.858 1.537 0.290 1.405 5.849
beta2_pelagic[12] 0.000 0.000 0.000 0.000 0.000
beta2_pelagic[13] 0.149 0.061 0.075 0.136 0.303
beta2_pelagic[14] 0.871 0.856 0.202 0.600 3.578
beta2_pelagic[15] 2.543 1.898 0.151 2.149 7.185
beta2_pelagic[16] 0.759 1.012 0.113 0.347 4.096
beta3_pelagic[1] 28.783 3.907 22.362 28.861 35.182
beta3_pelagic[2] 28.685 3.876 22.331 28.727 35.112
beta3_pelagic[3] 27.799 2.715 23.164 28.055 32.191
beta3_pelagic[4] 26.173 1.281 23.446 26.300 28.782
beta3_pelagic[5] 29.508 6.569 18.665 29.268 40.383
beta3_pelagic[6] 29.486 6.630 18.549 29.603 40.390
beta3_pelagic[7] 29.614 6.572 18.694 29.623 40.410
beta3_pelagic[8] 29.503 6.741 18.445 29.486 40.503
beta3_pelagic[9] 29.338 6.609 18.575 29.213 40.356
beta3_pelagic[10] 29.544 6.568 18.494 29.543 40.282
beta3_pelagic[11] 30.087 1.226 27.832 30.047 32.705
beta3_pelagic[12] 30.000 0.000 30.000 30.000 30.000
beta3_pelagic[13] 31.161 4.556 23.961 30.760 40.677
beta3_pelagic[14] 34.144 1.613 31.408 34.059 37.819
beta3_pelagic[15] 33.869 5.704 23.547 34.470 41.703
beta3_pelagic[16] 31.336 2.812 26.280 31.189 37.932
beta0_yellow[1] -0.483 0.080 -0.641 -0.483 -0.323
beta0_yellow[2] 0.545 0.080 0.379 0.551 0.682
beta0_yellow[3] -0.205 0.127 -0.449 -0.209 0.013
beta0_yellow[4] 0.884 0.257 0.123 0.952 1.186
beta0_yellow[5] 0.236 0.325 -0.448 0.261 0.793
beta0_yellow[6] 0.818 0.107 0.602 0.816 1.026
beta0_yellow[7] 1.014 0.100 0.821 1.014 1.209
beta0_yellow[8] 1.017 0.094 0.832 1.017 1.201
beta0_yellow[9] 0.443 0.099 0.252 0.442 0.638
beta0_yellow[10] 0.540 0.072 0.401 0.539 0.686
beta0_yellow[11] -0.200 0.127 -0.479 -0.188 0.014
beta0_yellow[12] -1.086 0.095 -1.267 -1.084 -0.900
beta0_yellow[13] -1.354 0.183 -1.778 -1.334 -1.059
beta0_yellow[14] -0.251 0.157 -0.566 -0.234 -0.008
beta0_yellow[15] -1.247 0.129 -1.546 -1.232 -1.030
beta0_yellow[16] -0.253 0.092 -0.434 -0.254 -0.071
beta1_yellow[1] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[2] 1.009 0.234 0.585 1.001 1.488
beta1_yellow[3] 0.461 0.247 0.017 0.475 0.924
beta1_yellow[4] 1.347 0.626 0.630 1.187 3.332
beta1_yellow[5] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[6] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[7] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[8] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[9] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[10] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[11] 0.302 0.368 0.001 0.121 1.208
beta1_yellow[12] 0.000 0.000 0.000 0.000 0.000
beta1_yellow[13] 1.307 1.297 0.028 0.851 5.181
beta1_yellow[14] 0.730 1.040 0.013 0.428 4.334
beta1_yellow[15] 0.164 0.196 0.011 0.085 0.670
beta1_yellow[16] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[1] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[2] -7.769 5.859 -21.918 -6.609 -0.470
beta2_yellow[3] -7.640 5.864 -22.076 -6.264 -0.440
beta2_yellow[4] -5.574 5.945 -21.034 -3.481 -0.116
beta2_yellow[5] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[6] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[7] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[8] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[9] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[10] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[11] -7.701 5.868 -21.568 -6.364 -0.472
beta2_yellow[12] 0.000 0.000 0.000 0.000 0.000
beta2_yellow[13] -5.826 5.324 -18.951 -4.026 -0.288
beta2_yellow[14] -7.215 5.604 -20.093 -5.943 -0.342
beta2_yellow[15] -7.666 5.582 -20.574 -6.576 -0.356
beta2_yellow[16] 0.000 0.000 0.000 0.000 0.000
beta3_yellow[1] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[2] 28.849 1.157 27.263 28.672 31.752
beta3_yellow[3] 32.349 3.279 23.378 32.556 38.386
beta3_yellow[4] 29.930 3.260 24.015 28.894 35.689
beta3_yellow[5] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[6] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[7] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[8] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[9] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[10] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[11] 34.926 3.680 30.333 33.852 43.292
beta3_yellow[12] 25.000 0.000 25.000 25.000 25.000
beta3_yellow[13] 33.323 2.797 29.196 33.367 39.911
beta3_yellow[14] 33.896 3.548 29.178 33.495 42.777
beta3_yellow[15] 36.372 3.745 30.323 36.177 43.444
beta3_yellow[16] 35.000 0.000 35.000 35.000 35.000
beta0_black[1] 0.012 0.077 -0.139 0.011 0.158
beta0_black[2] 1.855 0.031 1.792 1.855 1.916
beta0_black[3] 1.239 0.038 1.168 1.239 1.314
beta0_black[4] 2.426 0.048 2.333 2.427 2.520
beta0_black[5] 5.093 2.294 1.856 4.638 10.726
beta0_black[6] 5.265 2.206 2.635 4.689 11.054
beta0_black[7] 4.098 2.115 1.572 3.544 9.886
beta0_black[8] 0.897 0.220 0.505 0.887 1.347
beta0_black[9] 2.810 0.272 2.300 2.802 3.353
beta0_black[10] 1.287 0.050 1.193 1.287 1.384
beta0_black[11] 2.980 0.052 2.879 2.980 3.081
beta0_black[12] 4.909 0.192 4.534 4.907 5.290
beta0_black[13] 0.447 0.338 0.018 0.369 1.227
beta0_black[14] 2.855 0.106 2.645 2.857 3.059
beta0_black[15] 0.722 0.050 0.622 0.723 0.819
beta0_black[16] 3.771 0.072 3.630 3.771 3.915
beta2_black[1] 0.000 0.000 0.000 0.000 0.000
beta2_black[2] 0.000 0.000 0.000 0.000 0.000
beta2_black[3] 0.000 0.000 0.000 0.000 0.000
beta2_black[4] 0.000 0.000 0.000 0.000 0.000
beta2_black[5] 0.000 0.000 0.000 0.000 0.000
beta2_black[6] 0.000 0.000 0.000 0.000 0.000
beta2_black[7] 0.000 0.000 0.000 0.000 0.000
beta2_black[8] 0.000 0.000 0.000 0.000 0.000
beta2_black[9] 0.000 0.000 0.000 0.000 0.000
beta2_black[10] 0.000 0.000 0.000 0.000 0.000
beta2_black[11] 0.000 0.000 0.000 0.000 0.000
beta2_black[12] 0.000 0.000 0.000 0.000 0.000
beta2_black[13] -2.574 1.775 -7.034 -2.147 -0.317
beta2_black[14] 0.000 0.000 0.000 0.000 0.000
beta2_black[15] 0.000 0.000 0.000 0.000 0.000
beta2_black[16] 0.000 0.000 0.000 0.000 0.000
beta3_black[1] 25.000 0.000 25.000 25.000 25.000
beta3_black[2] 25.000 0.000 25.000 25.000 25.000
beta3_black[3] 25.000 0.000 25.000 25.000 25.000
beta3_black[4] 25.000 0.000 25.000 25.000 25.000
beta3_black[5] 25.000 0.000 25.000 25.000 25.000
beta3_black[6] 25.000 0.000 25.000 25.000 25.000
beta3_black[7] 25.000 0.000 25.000 25.000 25.000
beta3_black[8] 25.000 0.000 25.000 25.000 25.000
beta3_black[9] 25.000 0.000 25.000 25.000 25.000
beta3_black[10] 25.000 0.000 25.000 25.000 25.000
beta3_black[11] 25.000 0.000 25.000 25.000 25.000
beta3_black[12] 25.000 0.000 25.000 25.000 25.000
beta3_black[13] 37.343 3.416 26.611 38.144 40.981
beta3_black[14] 25.000 0.000 25.000 25.000 25.000
beta3_black[15] 25.000 0.000 25.000 25.000 25.000
beta3_black[16] 25.000 0.000 25.000 25.000 25.000
beta4_black[1] -0.326 0.091 -0.501 -0.325 -0.149
beta4_black[2] 0.258 0.053 0.155 0.258 0.362
beta4_black[3] -0.718 0.065 -0.846 -0.718 -0.593
beta4_black[4] 0.625 0.131 0.378 0.622 0.884
beta4_black[5] 0.818 1.823 -1.059 0.505 4.695
beta4_black[6] 0.864 1.651 -0.803 0.558 4.340
beta4_black[7] 0.814 1.561 -0.937 0.504 4.245
beta4_black[8] 0.078 0.290 -0.534 0.086 0.623
beta4_black[9] 1.024 0.812 -0.102 0.860 2.971
beta4_black[10] 0.073 0.065 -0.049 0.073 0.204
beta4_black[11] -0.609 0.067 -0.738 -0.609 -0.477
beta4_black[12] 1.236 0.290 0.697 1.226 1.815
beta4_black[13] -1.006 0.077 -1.158 -1.003 -0.862
beta4_black[14] 0.216 0.181 -0.142 0.216 0.577
beta4_black[15] 0.003 0.045 -0.084 0.003 0.091
beta4_black[16] -0.417 0.115 -0.646 -0.417 -0.193
beta0_dsr[11] -0.329 0.009 -0.347 -0.329 -0.310
beta0_dsr[12] 1.419 0.044 1.335 1.419 1.508
beta0_dsr[13] 0.155 0.021 0.112 0.154 0.194
beta0_dsr[14] -0.126 0.017 -0.160 -0.127 -0.092
beta0_dsr[15] 0.205 0.014 0.177 0.205 0.233
beta0_dsr[16] 0.161 0.016 0.129 0.161 0.193
beta1_dsr[11] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[12] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[13] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[14] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[15] 0.000 0.000 0.000 0.000 0.000
beta1_dsr[16] 0.000 0.000 0.000 0.000 0.000
beta4_dsr[11] 0.543 0.035 0.473 0.543 0.614
beta4_dsr[12] 0.799 0.226 0.374 0.796 1.279
beta4_dsr[13] -0.078 0.040 -0.157 -0.078 0.000
beta4_dsr[14] -0.209 0.072 -0.356 -0.209 -0.067
beta4_dsr[15] 0.217 0.022 0.172 0.217 0.260
beta4_dsr[16] -0.052 0.047 -0.145 -0.053 0.039
beta0_slope[11] -2.208 0.021 -2.250 -2.208 -2.168
beta0_slope[12] -5.328 0.212 -5.750 -5.325 -4.916
beta0_slope[13] -1.501 0.033 -1.565 -1.501 -1.437
beta0_slope[14] -2.711 0.051 -2.816 -2.710 -2.614
beta0_slope[15] -1.591 0.022 -1.634 -1.591 -1.550
beta0_slope[16] -2.651 0.041 -2.727 -2.653 -2.567
beta1_slope[11] 0.000 0.000 0.000 0.000 0.000
beta1_slope[12] 0.000 0.000 0.000 0.000 0.000
beta1_slope[13] 0.000 0.000 0.000 0.000 0.000
beta1_slope[14] 0.000 0.000 0.000 0.000 0.000
beta1_slope[15] 0.000 0.000 0.000 0.000 0.000
beta1_slope[16] 0.000 0.000 0.000 0.000 0.000
beta4_slope[11] -0.336 0.075 -0.494 -0.335 -0.192
beta4_slope[12] -4.680 0.803 -6.373 -4.636 -3.227
beta4_slope[13] 0.328 0.052 0.220 0.329 0.426
beta4_slope[14] -0.158 0.153 -0.469 -0.154 0.129
beta4_slope[15] -0.360 0.033 -0.427 -0.361 -0.295
beta4_slope[16] -0.141 0.104 -0.346 -0.142 0.058
beta5_slope[11] 3.064 0.105 2.869 3.062 3.283
beta5_slope[12] 1.999 0.384 1.316 1.980 2.819
beta5_slope[13] 0.881 0.081 0.722 0.879 1.043
beta5_slope[14] 3.743 0.361 3.075 3.734 4.491
beta5_slope[15] 1.308 0.069 1.178 1.306 1.445
beta5_slope[16] 2.727 0.239 2.281 2.721 3.221
tau_comp[1] 3.776 0.288 3.230 3.766 4.368
tau_comp[2] 2.855 0.261 2.376 2.844 3.389
tau_comp[3] 0.867 0.075 0.729 0.864 1.022
tau_comp[4] 1.780 0.277 1.282 1.761 2.386
tau_comp[5] 1.307 0.194 0.967 1.294 1.733
sigma_H[1] 0.209 0.050 0.119 0.206 0.316
sigma_H[2] 0.175 0.029 0.128 0.172 0.239
sigma_H[3] 0.187 0.040 0.116 0.184 0.276
sigma_H[4] 0.462 0.079 0.330 0.453 0.640
sigma_H[5] 0.965 0.203 0.596 0.951 1.399
sigma_H[6] 0.537 0.193 0.156 0.533 0.913
sigma_H[7] 0.309 0.062 0.215 0.300 0.452
sigma_H[8] 0.419 0.088 0.268 0.410 0.615
sigma_H[9] 0.512 0.119 0.329 0.497 0.788
sigma_H[10] 0.210 0.039 0.143 0.207 0.296
sigma_H[11] 0.266 0.044 0.193 0.262 0.364
sigma_H[12] 0.387 0.147 0.201 0.343 0.714
sigma_H[13] 0.218 0.037 0.153 0.215 0.297
sigma_H[14] 0.497 0.088 0.341 0.491 0.685
sigma_H[15] 0.231 0.037 0.167 0.227 0.310
sigma_H[16] 0.215 0.041 0.149 0.210 0.312
lambda_H[1] 2.254 3.392 0.078 1.075 11.748
lambda_H[2] 5.254 5.161 0.526 3.695 20.043
lambda_H[3] 5.233 7.616 0.264 2.709 26.319
lambda_H[4] 0.006 0.004 0.001 0.005 0.018
lambda_H[5] 2.929 7.362 0.025 0.688 23.772
lambda_H[6] 6.216 12.810 0.008 0.316 42.593
lambda_H[7] 0.012 0.008 0.002 0.010 0.032
lambda_H[8] 7.961 10.390 0.078 4.305 34.897
lambda_H[9] 0.015 0.010 0.003 0.013 0.040
lambda_H[10] 0.231 0.304 0.026 0.158 0.833
lambda_H[11] 0.224 0.317 0.010 0.105 1.097
lambda_H[12] 3.887 5.407 0.163 2.142 18.994
lambda_H[13] 3.222 3.057 0.263 2.342 11.118
lambda_H[14] 3.348 4.047 0.221 2.149 13.411
lambda_H[15] 0.019 0.020 0.003 0.014 0.063
lambda_H[16] 0.503 0.778 0.033 0.267 2.288
mu_lambda_H[1] 4.025 1.886 1.076 3.811 8.149
mu_lambda_H[2] 3.719 1.933 0.567 3.586 7.790
mu_lambda_H[3] 3.403 1.842 0.703 3.134 7.684
sigma_lambda_H[1] 8.037 4.412 1.735 7.209 18.182
sigma_lambda_H[2] 8.267 4.737 0.883 7.683 18.099
sigma_lambda_H[3] 6.275 4.039 0.960 5.417 16.315
beta_H[1,1] 6.504 1.313 3.400 6.709 8.402
beta_H[2,1] 9.763 0.610 8.442 9.798 10.851
beta_H[3,1] 7.832 0.780 5.995 7.938 9.123
beta_H[4,1] 10.639 7.707 -4.591 10.623 26.439
beta_H[5,1] -0.301 2.488 -5.716 -0.089 3.860
beta_H[6,1] 3.087 4.435 -7.738 4.697 8.375
beta_H[7,1] 1.365 5.511 -10.387 1.805 11.552
beta_H[8,1] 1.378 4.120 -2.623 1.232 3.652
beta_H[9,1] 13.668 5.772 2.340 13.617 25.583
beta_H[10,1] 6.807 1.859 2.845 6.894 10.336
beta_H[11,1] 4.849 3.552 -3.092 5.540 9.880
beta_H[12,1] 2.607 1.084 0.731 2.537 5.011
beta_H[13,1] 9.083 0.875 7.256 9.148 10.604
beta_H[14,1] 2.126 0.996 0.054 2.133 4.102
beta_H[15,1] -7.069 3.325 -13.330 -7.102 -0.162
beta_H[16,1] 3.879 2.977 -0.849 3.408 10.874
beta_H[1,2] 7.783 0.273 7.212 7.786 8.300
beta_H[2,2] 9.882 0.164 9.577 9.877 10.211
beta_H[3,2] 8.803 0.201 8.408 8.802 9.203
beta_H[4,2] 3.020 1.441 0.259 2.981 5.895
beta_H[5,2] 1.775 0.987 -0.154 1.788 3.640
beta_H[6,2] 6.020 1.090 3.592 6.187 7.795
beta_H[7,2] 2.391 1.027 0.502 2.318 4.534
beta_H[8,2] 2.989 1.155 1.317 3.139 4.274
beta_H[9,2] 3.175 1.095 1.088 3.138 5.354
beta_H[10,2] 8.051 0.366 7.298 8.062 8.770
beta_H[11,2] 9.839 0.634 8.897 9.731 11.223
beta_H[12,2] 3.935 0.369 3.277 3.922 4.719
beta_H[13,2] 9.156 0.242 8.703 9.149 9.666
beta_H[14,2] 4.022 0.352 3.342 4.018 4.731
beta_H[15,2] 11.555 0.604 10.281 11.573 12.708
beta_H[16,2] 4.451 0.810 2.969 4.422 6.113
beta_H[1,3] 8.442 0.264 7.966 8.435 8.992
beta_H[2,3] 9.855 0.145 9.560 9.854 10.140
beta_H[3,3] 9.466 0.160 9.156 9.462 9.803
beta_H[4,3] -2.169 0.865 -3.946 -2.138 -0.532
beta_H[5,3] 3.900 0.587 2.662 3.902 5.065
beta_H[6,3] 8.317 1.178 6.648 8.017 10.861
beta_H[7,3] -2.561 0.624 -3.828 -2.545 -1.392
beta_H[8,3] 5.308 0.549 4.710 5.214 6.432
beta_H[9,3] -2.732 0.713 -4.153 -2.710 -1.333
beta_H[10,3] 8.538 0.273 8.005 8.538 9.075
beta_H[11,3] 8.492 0.290 7.866 8.515 8.987
beta_H[12,3] 5.266 0.321 4.505 5.303 5.801
beta_H[13,3] 8.866 0.176 8.490 8.872 9.189
beta_H[14,3] 5.777 0.266 5.212 5.800 6.252
beta_H[15,3] 10.273 0.298 9.692 10.272 10.850
beta_H[16,3] 5.936 0.600 4.721 5.943 7.025
beta_H[1,4] 8.094 0.191 7.688 8.102 8.436
beta_H[2,4] 9.854 0.132 9.583 9.859 10.102
beta_H[3,4] 9.955 0.152 9.628 9.966 10.224
beta_H[4,4] 11.802 0.470 10.870 11.809 12.701
beta_H[5,4] 5.570 0.733 4.382 5.487 7.289
beta_H[6,4] 7.172 0.914 5.180 7.398 8.446
beta_H[7,4] 8.538 0.343 7.869 8.539 9.206
beta_H[8,4] 6.721 0.270 6.202 6.738 7.159
beta_H[9,4] 7.441 0.443 6.568 7.439 8.302
beta_H[10,4] 7.525 0.224 7.113 7.520 8.000
beta_H[11,4] 9.512 0.192 9.143 9.513 9.900
beta_H[12,4] 7.241 0.196 6.860 7.240 7.642
beta_H[13,4] 9.070 0.141 8.777 9.073 9.335
beta_H[14,4] 7.833 0.212 7.437 7.828 8.259
beta_H[15,4] 9.482 0.224 9.042 9.480 9.929
beta_H[16,4] 9.433 0.237 8.975 9.430 9.899
beta_H[1,5] 8.903 0.151 8.589 8.905 9.195
beta_H[2,5] 10.676 0.100 10.493 10.671 10.891
beta_H[3,5] 10.766 0.161 10.479 10.760 11.099
beta_H[4,5] 8.351 0.484 7.429 8.339 9.318
beta_H[5,5] 5.189 0.606 3.709 5.265 6.177
beta_H[6,5] 8.925 0.628 7.963 8.812 10.278
beta_H[7,5] 6.565 0.334 5.919 6.557 7.258
beta_H[8,5] 8.259 0.217 7.932 8.239 8.723
beta_H[9,5] 8.195 0.454 7.298 8.194 9.105
beta_H[10,5] 10.037 0.230 9.576 10.043 10.483
beta_H[11,5] 11.594 0.221 11.142 11.594 12.028
beta_H[12,5] 8.603 0.188 8.245 8.604 8.983
beta_H[13,5] 10.140 0.131 9.885 10.136 10.399
beta_H[14,5] 9.246 0.218 8.850 9.232 9.703
beta_H[15,5] 11.311 0.235 10.861 11.310 11.764
beta_H[16,5] 9.838 0.188 9.457 9.845 10.185
beta_H[1,6] 10.116 0.235 9.749 10.086 10.693
beta_H[2,6] 11.472 0.117 11.234 11.472 11.707
beta_H[3,6] 10.619 0.158 10.261 10.629 10.897
beta_H[4,6] 12.930 0.860 11.175 12.956 14.561
beta_H[5,6] 5.673 0.624 4.541 5.648 6.988
beta_H[6,6] 8.706 0.572 7.399 8.741 9.757
beta_H[7,6] 10.563 0.578 9.365 10.570 11.659
beta_H[8,6] 9.518 0.274 9.044 9.536 9.936
beta_H[9,6] 8.284 0.758 6.838 8.287 9.747
beta_H[10,6] 9.153 0.347 8.404 9.180 9.786
beta_H[11,6] 10.687 0.366 9.901 10.713 11.363
beta_H[12,6] 9.377 0.237 8.946 9.370 9.873
beta_H[13,6] 11.092 0.154 10.786 11.090 11.404
beta_H[14,6] 9.916 0.283 9.330 9.922 10.466
beta_H[15,6] 10.552 0.414 9.743 10.538 11.354
beta_H[16,6] 10.606 0.254 10.081 10.617 11.094
beta_H[1,7] 10.194 1.224 7.037 10.470 11.869
beta_H[2,7] 12.113 0.523 10.986 12.149 13.057
beta_H[3,7] 10.299 0.681 8.822 10.349 11.478
beta_H[4,7] 1.340 4.448 -7.221 1.234 10.295
beta_H[5,7] 5.858 2.051 1.749 5.878 10.100
beta_H[6,7] 7.101 3.189 -0.764 8.282 10.878
beta_H[7,7] 7.575 2.955 1.937 7.507 13.595
beta_H[8,7] 10.885 0.990 9.237 10.870 12.531
beta_H[9,7] 5.010 3.895 -2.732 5.150 12.677
beta_H[10,7] 10.008 1.624 7.153 9.836 13.721
beta_H[11,7] 11.400 1.816 8.095 11.214 15.363
beta_H[12,7] 9.795 0.979 7.629 9.867 11.511
beta_H[13,7] 11.726 0.702 10.102 11.780 12.912
beta_H[14,7] 10.653 0.925 8.736 10.696 12.345
beta_H[15,7] 13.978 2.189 9.755 14.026 18.237
beta_H[16,7] 12.807 1.369 10.461 12.693 15.754
beta0_H[1] 8.891 17.591 -26.689 8.896 44.445
beta0_H[2] 10.502 7.972 -5.421 10.517 25.952
beta0_H[3] 9.440 10.357 -11.963 9.629 30.120
beta0_H[4] 12.337 183.186 -348.316 12.262 373.611
beta0_H[5] 3.636 27.006 -54.516 4.405 56.307
beta0_H[6] 9.283 51.206 -102.918 7.969 121.501
beta0_H[7] 5.508 142.891 -290.794 4.494 290.885
beta0_H[8] 6.328 34.385 -19.249 6.355 32.966
beta0_H[9] 4.805 121.092 -230.356 3.090 253.128
beta0_H[10] 7.652 37.289 -73.394 8.284 82.453
beta0_H[11] 10.256 52.911 -98.632 9.977 122.540
beta0_H[12] 6.748 12.563 -17.128 6.849 30.952
beta0_H[13] 9.731 10.532 -9.761 9.660 30.742
beta0_H[14] 6.832 11.546 -15.605 6.926 30.000
beta0_H[15] 4.877 116.679 -244.180 7.245 236.494
beta0_H[16] 9.029 31.187 -51.965 8.566 71.765